Viscoelastic fluids in thin domains: a mathematical proof 
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Abstract 



The present paper deals with non Newtonian viscoelastic flows of Oldroyd-B type in thin 
domains. Such geometries arise for example in the context of lubrication. More precisely, we 
justify rigorously the asymptotic model obtained heuristically by proving the mathematical 
convergence of the Navier-Stokes/ Oldroyd-B sytem towards the asymptotic model. 
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1 Introduction 

m ■ 

■ This paper concerns the study of a viscoelastic fluid flow in a thin gap, the motion of which is 
imposed due to non homogeneous boundary conditions. 

When a Newtonian flow is contained between two close given surfaces in relative motion, it is 
well known that it is possible to replace the Stokes or Navier-Stokes equations governing the fluid's 
motion by a simpler asymptotic model. The asymptotic pressure is proved to be independent 
^ of the normal direction to the close surfaces and obeys the Reynolds thin fllm equation whose 

^ ' coefficients include the velocities, the geometrical description of the surrounding surfaces and 
some rheological characteristics of the fiuid. As a following step, the computation of this pressure 
allows an asymptotic velocity of the fiuid to be easily computed. Such asymptotic procedure 
first proposed in a formal way by Reynolds [2] has been rigorously confirmed for Newtonian 
stationary fiow [1], and then generalized in a lot of situations covering numerous applications for 
both compressible fiuid [14], unsteady cases [3], multifiuid fiows [15]. 

It is well known however that in numerous applications, the fiuid to be considered is a non 
Newtonian one. This is the case for numerous biological fiuids, modern lubricants in engineering 
applications due to the additives they contain, polymers in injection or molding process. In all of 
these applications, there are situations in which the fiow is anisotropic. It is usual to take account 
of this geometrical effect in order to simplify the three-dimensional equations of the motion. 
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trying to recover two dimensional Reynolds like equation with respect to the pressure only. Such 
procedures are more often heuristic ones. Nevertheless, some mathematical works appeared in the 
literature to justify them. They include thin film asymptotic studies of Bingham flow [9], quasi 
Newtonian flow (Carreau's law, power law or Williamson's law, in which various stress-velocity 
relations are chosen: [7], [6], [16]) and also micro polar ones [5]. It has been possible to obtain 
rigorously some thin film approximation for such fluids using a so called generalized Reynolds 
equation for the pressure. 

However in the preceding examples, elasticity effects are neglected. Introduction of such vis- 
coelastic behavior is characterized by the Deborah number which is related to the relaxation time. 
One of the most popular laws is the Oldroyd-B model whose constitutive equation is an interpola- 
tion between purely viscous and purely elastic models, thus introducing an additional parameter 
which describes the relative proportion of both behaviors. A formal procedure has been proposed 
in [4]. However, the asymptotic system so obtained lacks the usual characteristic of classical 
generalized Reynolds equation as it has not been possible to gain an equation in the asymptotic 
pressure only. Both velocity u* and pressure p* are coupled by a non linear system. 

It is the goal of this paper to justify rigorously this asymptotic system. Section 2 is devoted 
to the precise statement of the 3-D problem. One difficulty has been to find an existence theorem 
for the general Oldroyd-B model, acting as a starting point for the mathematical procedure. Most 
of the existence theorems, however, deal with small data or small time assumptions. To control 
this kind of property with respect to the smallness of the gap appears somewhat difficult. So we 
are led to consider a more particular Oldroyd-B model, for which unconditional existence theorem 
has been proved [13]. Moreover, a specific attention is devoted to the boundary conditions to be 
introduced both on the velocity and on the stress. The goal is to use "well prepared" boundary 
conditions so as to prevent boundary layer on the lateral side of the domain. 

In Section 3, after suitable scaling procedure, asymptotic expansions of both pressure, viscosity 
and stress are introduced, taking into account the previous formal results from [4]. Section 4 is 
mainly concerned with the proof of some additional regularity properties for the formal asymptotic 
solution. Assuming some restrictions on the rheological parameters, it will be proved that it is 
possible to gain a C'^ regularity for p* , k > 1, which in turn improves the regularity of u* and the 
stress tensor a*. This result is obtained by introducing a differential Cauchy system satisfied by 
the derivative of p* . Finally, section 5, is devoted to the convergence towards zero of the second 
term of the asymptotic expansions, which in turn proves the convergence of the solution of the 
real 3-D problem towards u*, p*, a* (Theorems 5.4 and 5.6). 

2 Introduction of the problem and known results 
2.1 Formulation of the problem 

We consider unsteady incompressible flows of viscoelastic fluids, which are ruled by Oldroyd's law, 
in a thin domain = {{x,y) G M",a; G co and < y < £h{x)}, where lo is an (n— l)-dimensional 
domain, with n = 2 or n = S {x = xi or x = {xi,X2)), as in Figure 1. 
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Figure 1: Domain fig 

The following hypotheses on h are required: 

Vx G Col, < /lo < < and /j'^gC^(cj). 

Let u'^ = ("uf , ul, "uD be the velocity field in the three-dimensional case, or = (ufjul) in the 
two-dimensional case, the pressure, and <t^ the stress symmetric tensor in the domain Cl'^ . Bold 
letters stand for vectorial or tensorial functions, the notation / corresponds to a function / defined 
in the domain (l^, and the superscript ^ denotes the dependence on e. 

Formulation of the problem The following formulation of the problem holds in (0, cxd) x Cl^: 

p dfU^ + p u'^ ■ V-u^ — (1 — ^)'^ ^"i*^ + = V • <T^ , 

< V-u^ = 0, (2.1) 

A {dta" + ■ Vo-^ + g{&'',Vu'')) + a" = 2ruD{u^) , 

where the nonlinear terms g{a-^ jVii^), the vorticity tensor W{u^) and the deformation tensor 
D{u^) are given by: 

g{a^, Vu^) = -W{u^) ■ + ■ W{u^), 
W{u') = and D{u') = . 

In this formulation, the physical parameters are the viscosity u, the density p, and the relaxation 
time A. The parameter A is related to the viscoelastic behavior and the Deborah number. The 
parameter r G [0, 1) describes the relative proportion of the viscous and elastic behavior. 
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Initial conditions This problem is considered with the following initial conditions: 

u^\t=o = ul, &^\t=o = o-g, (2.2) 

for Uq G L'^{Ci^), <Tq G L^{Cl'^). The bold notation L'^{Cl^) denotes the set of vectorial or tensorial 
functions whose all components belong to L^(f2^). 

Boundary conditions Dirichlet boundary conditions are set on top and bottom of the domain, 
and the conditions on the lateral part of the boundary r|^, defined by 

tl = {{x,y) G M",a; G 5a; and < y < £h{x)} , 

will be specified later (in section 4.2). Therefore, it is possible to write the boundary conditions 
in a shortened way: 

^lane = (2-3) 

where is a given function such that G H^/'^{d()F) and satisfying J'^|j^=/je = 0, J'^|y=o = (s, 0). 
This function will be fully determined in Subsection 4.2. 

Since <t^ satisfies a transport equation in the domain Cl^ , it remains to impose boundary conditions 
on (T^ on the part of the boundary where tt^ is an incoming velocity. Let us define the part of 
tl such that >|pe^ • n < 0, and f i = ff, \ f ^. We set 

«T"|f^=^", (2.4) 

where 6^ is a given function in i3"^/^(r^) which will also be determined in Subsection 4.2. 
Moreover, since the pressure is defined up to a constant, the mean pressure is chosen to be zero: 
j' f = Q. 

Notations Let us introduce the following function space: 

v = [<peH^{n'),v-cp = o'\, 

and the following notations, that will be used in the following. For / defined in 0,^: 

• 1/1 denotes the L^-norm in Cl^, 

• \f\p denotes the L^'-norm in J^^, for 2 < p < +oo, 

• the spaces C"^('tl^) for m > 1 are equipped with the norms H/Hc" = |/|oo + Yl |/^*^|oo- 

i=l 

For / defined in R"^ x D,^ , ||/||l«(l/3) denotes the norm of the space L"(0,oo,L^(fi^)), with 1 < 
a, P < oo. 
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2.2 Existence theorem in the domain Cl^ 



Theorem 2.1. Fore > fixed, problem (2.1)-(2.3) admits a weak solution 



Proof. This result is proved in [13]. 



□ 



Remark 2.2. Let us emphasize that for the following, it is essential to know the global (in time) 
existence of a solution for problem (2.1)-(2.3). Other existence theorems have been proved for this 
problem, for example in [12], [11], ]10], but these theorems are either local in time (on a time 
interval [0, T^]), or a small data assumption is needed. In this work, these theorems cannot be 
used, since there is no control on the behavior ofT^ (or equivalently of the data) when e tends to 
zero, in particular may tend to zero. 

Consequently, this work is restricted to the specific case treated in [13], taking one parameter of 
the Oldroyd model to be zero. In all generality, the non-linear term reads g{a,'Vu) = —W{u) ■ 
a + cr ■ W{u) — a{cr ■ D{u) + D{u) ■ cr), which is called objective derivative. Here the parameter 
a is taken to be zero. This case corresponds to the so-called Jaumann derivative. 

Remark 2.3. The following computations are made in the two-dimensional case (i.e. uj = (0, L) is 
a one-dimensional domain) for the sake of simplicity. However, note that except for the regularity 
obtained for the limit problem in Section 4-3, all estimates are independent of the dimension, thus 
the corresponding computations should apply to the three-dimensional case. 

Regularizing the system In the proof of the preceding theorem, the existence of a solution is 
achieved by regularization. Therefore, this study only concerns solutions obtained as the limit of 
a regularized problem approximating (2.1), in which an additional term —rjA&^^ is added to the 
Oldroyd equation, with 77 > a small parameter. Here a regularization of the form —r]A(^&^'^ — G) 
is chosen, with G a symmetric tensor in H'^{i).^) independent of rj and e which will be precised 
later. After obtaining the needed energy estimates uniformly in r], we will let r] tend to zero. This 
approach allows to multiply the Oldroyd equation by a^'^', since a-^"^ is regular enough. Of course, 
one can choose another regularization which leads to energy estimates which are uniform in the 
regularization parameter. 

Furthermore, because of the regularizing term, boundary conditions on the whole boundary are 
needed. Let us write o-^^lg^^ = 6^'^, where 0^'^ is now a function of H^/'^(dCl^), which will be 
determined later by equation (4.3). 
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3 Asymptotic expansions 



3.1 Renormalization of the domain 

y 

Introducing a new variable z = the system (2.1) can be rewritten in a fixed re-scaled domain: 

ri = {(x, z) £ M", X £ui and < z < h{x)} . 

For a function / defined in / is defined in Q by f{x,z) = f{x,ez). For a function / G 

\f\p still denotes the I^-norm in J7, and similar notations hold for the other norms. Moreover, the 

regularizing term rjAcr^'^ is introduced. Denoting a^^ = | er/ I ' ^^^'^ similar notations for 



the components of G, it holds in (0, oo) x Q, : 

p Stul"^ -{l-r)u A,nf ' + d^p'^ - d^a^i - ^d^al^^ = , 

p 5tu^> -{l-r)v A,ul'> + U,p'^ - d,a{l - U.all = , 
• u^'? = , 

A 



Stall - N{u'\ al'D) + a\l - r]A,{alJ - Gu) - 2ri^d,uT = , 

(3.1) 

where the convective derivative 6t is given by St = df + u^^ ■ V^. The derivation operators are 
defined as follows: = ^dx, = H — 2^2" non-linear terms N are given by 



^{'f^^ /) = ( 9xU2 - ^d^ui ) /. 
3.2 Asymptotic expansions 

It has been proposed in [4] that when r), e tend zero, {u^'^ , p^'^ , cr^'^) tends formally to a triplet 
{u* ,p* , a*) satisfying a system that will be given later in (4.1). This analysis leads to the intro- 
duction of the following asymptotic expansions: 

= ul + vY' and = + , (3.2) 

rr^r^ = Ig-* + l,-^^, (3.4) 



£ £ 



With a* = [ ""1} '^f ) , and r^^' = [ ^fj . If denoting u* = (^xj,u*), and v^^ = {vY,vT), 

\^12 ^22 J \''l2 "^22 / 

(3.2) becomes u^"^ = u* + v^'^ . 

The scaling orders chosen for the pressure and the different components of the velocity field and 
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of the stress tensor are motivated by some mathematical and physical remarks. Classically, the 
pressure has to be of order if the horizontal velocity is of order 1 (see [2] for the rigorous 
explanation). On the other hand, the stress tensor has to be of order 1/e and the Deborah number 
A of order e in order to balance the Newtonian and non-Newtonian contribution in Oldroyd 
equation (see [4]). Hence; let A = eA*. 

A wise choice of the function G in the regularizing term is G = <t*. The regularity of G in 
-ff2(0) is proved by Theorem 4.4 (where it is proved that dla* G C°(il), d,^dz(T* £ C°(f7) and 
d'^a* e e^{n), thus Aa* e 2.2(0)). A formal substitution of (3.2), (3.3), (3.4) in (3.1) leads to 
the following system: 

p dtvl^ - (1 - r)^ A.^" + \d,q'^ - -d^rl^ - ^d.rl^ = + -C, + \c[, 

p dtvl'' - (1 - r)uA,v'2' + - ^d^rt^ - ^d^r^^ = ^^2' + + ^C^, 

V • v^" = V • u*, 

A* (d,ri7 - N{v^^ rt^)) + ^r^^ - r^A.r^^ - 2r,.d,vr' = L^^ + 

(3.5) 

with the following notations: dt = dt + v^'^ ■ V is the so-called convective derivative, the non-linear 

terms N{v^^,f) = ^dxV2^ dzvY'^ / for / G L^(p,) and the following linear (with respect to 

v^'^) and constant terms 

L^^ = -p -^ul-pu* ■ Vv^ -p dtu\ - p w* • Vtxt + (1 - r)vdlu\, 
C{ = {\- T)vdlu{ - d^p* + d,al2\ 



- pe dtu*2 -pen*- Vu^ + £ (1 - r)vdi,u*2 + (1 - r)vd^u*2 + 5j;(7*2, 
C2 = dzcr22, 
C'2 = d,p*. 

For the Oldroyd equation, the following linear (with respect to v and r) and constant terms 
appear: 

Lfl =Cl1 + A* {-dtal^ - u* ■ Vfj^i + ed^u*2al2) + 2riyd^ul, 
with Cfl = A* {ed^ulrll + e^a^uf (7^2 - t^^'' • VaJi - u* ■ VtID , 
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with = -A* {ed^u*2Tf^ + eda,V2al2 + v'" ■ Va^s + • Vt^^) , 



'-22 



■'-12 

A* 

- y (25tcri2 + 2w* • V(7^2 + 5x^2 (f^ii - 0-22)) + rv£dxu*2, 
~ A* A* 

Lfs" = - Y ((^.r"I(^n - 7-I2 ) + - ^22)) + Y^z^^iKi - 0-22) - <^\2 + r-i^S^tXi; 

^ V ' 

•"-12 

Note that the first order derivatives of <t* occur in the terms L^"^ and C^^. It will be shown in 
Theorem 4.4 that cr* has sufficient regularity. 

Let us observe also that equations (3.5) are similar to (3.1), except for the linear terms on the 
right. Thus the energy estimates will be obtained similarly for both systems, multiplying Navier- 
Stokes equation by the velocity and Oldroyd equation by the stress tensor, and integrating over 



4 Limit equations 
4.1 Limit system 

In an heuristic way, the following system of equations satisfied by w*, p*, (T* is infered from (3.5): 
It*, p*, <T* are steady functions solutions of: 

' (1 - T)vdlu\ - d^f + d,Gi2 = 0, 

^^V* = 0, 

V • w* = 0, 



\*d,ulal2 + (tI, = 0, 
A* 

— ^dzulial^ - 0-22) + = rud^ul, 
-\*dzulal2 + 0-22 = 0. 



(4.1) 
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This system is equipped with the following boundary condition (Dirichlet condition on the upper 
and lower part of the boundary, flux imposed on the lateral part of the boundary): 



tt* = , for ^; = h{x), 

u* = {s,0), for2; = 0, ^4 2) 

h{x) 

J u*dz-n = ^o on r^. 
>^ 

The compatibility condition reads / $o = 0- Moreover, since p* is defined up to a constant, the 

du) 

mean pressure is taken to be zero: J p* = 0. 

a, 

Remark 4.1. Each equation of the preceding system (4.1) is obtained by cancelling the constant 
part (i.e. the part independent of v^'^ , q^'^ , t^'^) of respectively C[, C'2, V • u* , L'^^ , L'l2 , L22. 

4.2 Determination of the boundary conditions 

Remark 4.2. The lateral boundary conditions on u* do not depend on the ones on u^^, but only 
on the flux. Therefore, different boundary conditions on u^'^ corresponding to the same fl,ux lead to 
the same limit problem. This is a classical fact when passing from a two-dimensional problem to 
a one- dimensional problem (or similarly from a three-dimensional problem to a two-dimensional 
one), and has already been observed in [2] for example. Here, in order to avoid boundary layers, 
u^^ = u* is imposed on the lateral part of the boundary. 

Similarly, any value of cr'^'^ on the boundary leads to the same limit problem. Again, in order to 
avoid boundary layers, well-prepared boundary conditions are also chosen for cr^'^ . 

The preceding remark allows to define precisely the function introduced in (2.3). Since 
u*\r^ e i?^/^(rL), it is possible to construct G H^^'^{d^l) satisfying J'^\z=h = 0, J^\z=o = (s,0) 
and J^Ifi^ = u*\rj^. Therefore, the boundary conditions on u^''' become 

u^'^ = , for z = h{x), 
u^v = (s^ 0) , for z = 0, 
u^'' = u* on r^. 

Thus u^^Iqq = u*\qq, and v^'^ will satisfy zero boundary conditions: v^'^\qq = 0. 

Moreover, since cr* G H^(Q,) (see Theorem 4.4 for this regularity result), 0^ can be defined as 

follows: 

6>^ = tT*|r+ eif^/2(r+). (4.3) 

Therefore 

cr^''|r+=o-*|r+, 

and this implies that r'^'^lr , = 0. 
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On the other part r_ of the boundary, cr^'^ is chosen such that tr^^ • n|r_ = f * • ?^|r_ , for example 
tr^'J |r_ = (T*\r_. 

4.3 Existence of a solution to the limit problem 

System (4.1)-(4.2) has aheady been studied in [4]. 

Theorem 4.3. Assume that r < 8/9. Then system (4.1)-(4.2) has a unique solution satisfying 
u* e L'^in), d^u* e p* e H^{u), a* e (4.4) 

Proof. This result has been proved in [4]. □ 

This existence result is not sufficient for this study. Therefore, the following stronger regularity 
result is proved on the limit problem (4.1)-(4.2). 

Theorem 4.4. Assume r < 2/9. If h e H^{u), for k G N*, then the unique solution {u*,p*,cr*) 
of the system (4.1)-(4.2) satisfies 

p* G C''+\id), ul, d,ul, dlul G C^+^{9), a*, d,a* G e''+\n), 

d^ul G C''{n), u*2, d,ul,dlul G C^(f^), 9^0-* G C^ip), (4.5) 
d^ul G C*^-^^^)- 

Proof. Let us observe that system (4.1) can be expressed as a system on u\, p* only. Using (4.1), 

a\-y, (722 be expressed as functions of cr*2 and dzu\. Indeed, from the fourth and the last 
equations of (4.1), it holds that 

o"22 = -0-11 = OzU'^ai'i- (4.6) 

Moreover, the divergence-free equation can be rewritten in order to eliminate Ug. Integrating this 
equation between z = and z = h, and using the fact that «2U=o = m2=/i = u\\z=h = 0, it 
follows: 

d^{^j^u\d^=Q. (4.7) 
Thus, the system in u\, p* can be written in the following form: 



12 1 + A*2|a^^x*|2' 



— v{l — r)d1u\ — dz<7\2 + (^xP* = 0, with (t\2 
dzP* = 0, 

h 



(4.8) 



djjuldz\=0, 



equipped with the boundary conditions stated in (4.2) and the condition J p* = 0. 

Q 
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For the sake of readability, the superscripts * are omitted in the rest of this section. 

Denote q = dxP- Let G C°°(M) defined by (f){t) = v{l — r)t + ry^ ■ The first equation of 

(4.8) becomes q = dz{(l){dzU\)). 

A simple study of function ^ allows to show the following properties: 

< ( 1 - ^ ) < l^'loo < z^, and Mt) > ±oo. (4.9) 

y 8 y t— »±oo 

Therefore the function is invertible, and ijj = belongs to C°°(M). Moreover, V' is an increasing 
function as (p. Integrating q = dzi<p{dzUi)) with respect to z between and z, the first equation 
of (4.8) becomes: 

(j){dzUi{x, z)) = q{x) z + k{x), 
where k{x) is a integration constant. Therefore, it follows that 

dzUi{x, z) = il^iqix) z + k{x)). 

Since "^11^=0 = s, the integration between and z of the preceding equation yields: 

ui{x,z)=s+[ 'ilj{q{x)t + K{x))dt. (4.10) 
Jo 

The boundary condition tti|/i(a;) = implies also: 

[ \{q{x)t + k{x)) + s = 0. (4.11) 

rh 

For {h, q, s, k) G M^, let us introduce F{h, q, s,k) = / il^{qt + k) + s. 

Jo 

Lemma 4.5. For any {h, q, s) G there exists an unique k G R such that F{h, q, s, k) = 0. 



Proof. • If such an k exists, it is unique from the implicit function theorem, since for all 
{h,q,s,K) G R'^ 

— {h,q,s,K)= / tp'{qt + Kjdt > 0. 
oi^ Jo 



• The following limits are computed, using the fact that lim ip(t) = ±00: 

t— >±oo 

lim F{h, q, s, k) = +00 and lim F{h, q, s, k) = —00. 

K—I- + 00 K—^ — CO 

Therefore, there exists k G M such that F{h, q, s, k) = 0. Let us denote K{h, q, s) = k. By 
the imphcit function theorem, K G C°°(M^). 

□ 
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Therefore, the following expression holds for {h,q,s) G M^: 

F{h,q,s,K{h,q,s))=0. (4.12) 

It is now possible to obtain an information on the sign of dqK. Indeed, deriving the expression 
(4.12) with respect to q, it follows 

dgF + d^F dgK = 0. 

ph rh 

For ^ > 0, since d^F = / til)'{qt + K)dt > and daF = / ^'(gi + K)dt > 0, dgK is strictly 

Jo Jo 

negative. 

Now, using equation (4.7) and the expression (4.10) for u, it follows: 

J^^^ da,(i^{q{x)t + K{h{x),q{x),s))ytdz = 0. 
or if changing the direction of integration 

^ ^ \h{x) - t)d^(^^{q{x)t + K{h{x),q{x), s))'^dt = 0. 
This can be rewritten as 

q'{x) ^ ^ {h{x) - t) {{t + dqK{h{x), q{x), s)) ^j' {q{x)t + K{h{x), q{x), s)) dt 

= - f^^'\Kx) - t) (h'{x)dhK{h{x),q{x), s)) ^'{q{x)t + K{h{x), q{x), s)) dt, 
which can be seen as an ordinary differential equation in q. Let 

U{x,q)= ^ ^ (h{x)-t^ (^t + dqK{h{x),q,s)^i^'{qt + K{h{x),q,s))dt, 
V{x,q)= ^ ^ {h{x)-t) (h'{x)dhK{h{x),q,s)^ i;' {qt + K{hix),q, s))dt. 

The differential equation becomes U {x , q{x)) q' {x) = —V{x,q(x)) for x E co. Note that this 
equation is in some sense a generalized Reynolds equation for the pressure. 

Lemma 4.6. Let r < 2/9. Then U{x, g) < for any (x, q) e lo xR. 

Proof. Let {x,q) € a; x R. Equation (4.11) and the definition (4.12) of K imply: 

rh{x) 



[ tp{qt + K{h{x),q,s))dt = -s, 
Jo 
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which becomes, after derivation with respect to q 

fh(x) 



(t + dgK{h{x),q,s))^'{qt + K{h{x),q,s))dt = 0. (4.13) 



With the notation K'(x,q) = dqK(h{x),q, s), (4.13) imphes 

rh{x) 



I t'ilj'{qt + K{h{x),q,s))dt 

/ ip'{qt + K{h{x),q,s))dt 
Jo 



Now, using this expression, U{x,q) can be simphfied 

fhix) 



U{x,q)=j^ -t(t + dqK{h{x),q,s))ilj'{qt + K{h{x),q,s))dt. (4.14) 



Recalhng the estimate of l^loo in (4.9), it follows that for any t G R: 



V </>'(V'(i) iy{l-9r/8) 



Let m = -, M = — — ^—rrr. Then 
u' 1/(1 -9r/8) 



^^^""^ <K'ix,q)< 



2m - ^ - 2M 

Now, (4.14) implies that: 

^ m M 
In order to prove that U remains strictly negative, it suinces to prove that < i.e. that 

m 3 2 

— > -, which is satisfied under the condition r < -. □ 
M 4 9 

It is possible to apply Picard-Lindelof theorem (or Cauchy-Lipschitz theorem) to the ordinary 
differential equation —U{x,q{x))q'{x) = V{x,q{x)), as U remains strictly negative by Lemma 
4.6. Since and K are C°°-functions, the regularity of q' is determined by the regularity of q and 
h. By hypothesis, h belongs to H^ioj), with A; G N, hence h G LP'{oj). Moreover, Theorem 4.3 
implies that q G i^^(w). Thus q' G L?{i^), which means q G H^{uj). 

Iterating this process as long as h is regular, h G H^{uj) and q G H^{oj) implies that q' G H^{oj), 
thus dxP = q & H'^^^{io), and p G H^^'^{u). By the classical Sobolev embedding, p belongs to 
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Last, recalling the expression (4.10), it follows that ui G C'^+-'^(a)), and, taking the first and second 
derivatives of (4.10) with respect to z, that dzUi, d^ui also belong to C^'^^{ijj). 
As observed in the introduction of the proof, <t and U2 are given as functions of p, ui, and the 
needed regularity follows. □ 

Remark 4.7. Since in practical applications, h is very regular (h G C°°{Cd) ), the preceding theorem 
gives as much regularity as wanted. In particular, the following result will he useful subsequently. 

Corollary 4.8. Assmne r < 2/9. If h & II^{(jj), then the unique solution {u*,p*,(T*) of the 
system (4.1)-(4.2) satisfies 

d:,u\ G C^ip), u*2, d^u*2,dlu*2 e C^ip), d^^a* G C\n), (4.15) 
Proof. It suffices to take A; = 1 in the preceding theorem 4.4. □ 

5 Convergence of the remainders 
5.1 Equations on the remainders 

From now on, the superscript '^'^ are dropped although the functions still depend on e and rj. Using 
the equations (4.1), system (3.5) becomes 

p dtvi - (1 - r)uA^vi + \d:^q - -d^Tn - = Li + -Ci, (5.1a) 

1 \ 1 1^1 
p dtV2 - (1 - r)i/AeV2 + -jd^q 2 ^^"^12 ^dzT22 = + -3C2, (5.1b) 

V-v = 0, (5.1c) 
A*dtrii - \*N{v,Ti2) + -Til - V^eTu - 2rud^vi = Lu + -L'^^ + 7]A,a*n, (5. Id) 

X* 1 / 1 \ 1 

A*dtri2 + Y^^'"' '^^^ ~ '^22) + -7-12 - V^eTi2 - ri^ I dxV2 + -9zVi j = L12 + -L'y2 + ??^£0-i^>5.1e) 

1 2rv 1 
A*dtr22 + A*iV(v, ri2) + -T22 - r]A^T22 ^"^^^2 = 1^22 + --^^22 + n^ecr22, (5-lf) 

with the new quantities 

Li ^ Ci ~ p u* ■ V?ii + (1 — r)vd1u\, 

L2 = C2- p £^u* ■ Vul + (1 - r)vdlul + (1 - r)z/9^n^ + d^al^, 
Lii = Cii + A* {-u* ■ Vail + £dxU2crl2) + 2rudxUi, 

T' — r' 

-^11 — ■'-11) 

A* 

^12 = >Ci2 - Y (2tt* • Vcr^2 + dxU*2{ali - (722)) + rved^ul, 

T' — r' 

-^12 — ■'-125 
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= >C22 - A* {u* ■ Vc722 + ed^U2al2) + 2^8 zul, 
T' — r' 

and with the initial and boundary conditions 

v\t=o = UQ-u*, T\t=o = (JQ - a* , = 0, T|r+ = 0. (5.2) 

Let us observe that both initial conditions ^14=0 and T|t=o belong to L^(J7). v, q and r are defined 
by (3.2), (3.3), (3.4). Prom the existence theorem 2.1 for (u,p, tr) and theorem 4.3 for {u*,p*,(T*), 
it follows that system (5.1) admits a solution {v,q,T) G L'^{0,oo, H^{^)) x L^(0, oo, L^(0)) x 
C{0,oo,L'^{n)) for r < 8/9. 

5.2 Convergence of v and r 

Before starting the a priori estimates, let us explain how the non-linear terms in (5.1) are han- 
dled. The non-linear terms v ■ Vv of Navier-Stokes equation and v ■ Vr of Oldroyd equation 
are treated with the following Lemma 5.1. On the other hand, the non-linear terms N(v,t) = 
{sd^V2 - Id.vi) T in (5.1d)-(5.1f) are zero when multiplied by r. 

Lemma 5.1. Let n be the exterior normal of the domain O,. Let (j) G H^{0,) be a vector field 
satisfying V • = and (p ■ n\QQ = 0. Let w G H^{Vt). Then 



(p ■ Vw w = Q. 



Proof. By integration by parts: 



(j) ■ Vw w = — J V • ■w'^ — J (p ■ Vw w + J (f) ■ n = 0. 

a =0 n an =o 



The classical approach consists in obtaining a priori estimates for v. 



□ 



Proposition 5.2. Let {v,q,T) be a solution of (5.1). Then v = {vi, V2) satisfy the following 
inequality for e small enough: 

d 3 

rvp- + \ev2\'') + -r(l - (|Ve^;i|' + \eVeV2\^) < -Vi -V2 + C, (5.3) 

2ru f 2rv f f 2rv f 

where T>i = j rn dxVi + \ t\2 Qz^x, V2 = 2rv / tu dxV2 + / r22 dzV2 and C is a 

constant independent of e. 

Proof. The proof consists in obtaining classical a priori estimates on both vi and V2. 
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Step 1. Let us multiply (5.1a) by vi and integrate over Observe that vi is regular enough 

to do so. Since v\gQ = 0, the boundary terms in the integration by parts are all zero. For example 

— J AgViVi = JlVgi^ip. Moreover, the convection terms / v ■ VviVi contained in / dtvivi are 

ci n n n 

equal to zero by Lemma 5.1, since V • w = and v\gci = 0. It follows: 

n a 

(5.4) 



-Vi/2ru 



It remains to estimate the terms / Livi and / CiVi. 

n n 

Main idea Estimates of the form: J Livi + ^ f Ci Ui < C+KilV^fi P + 1 92^2 P will be proved, 

n h 

where C is a constant independent of e and where the constants ki, K2 satisfy ki,K2 < (1 — r)z^/4. 
These constants will be precised later in the proof. 

In the following, C, q and Mj will denote some constants independent of e and rj, which might 
depend on on the physical parameters of the problem and on u* , cr* in sufficiently regular 
norms. 

• Let us estimate first the linear (with respect to v) term jCi of Li. To this end, Poincare inequality 
is useful: for / G L^(J7), with f\z=h = 0, |/| ^ Cp\dzf\. The constant Cp only depends on Q. 



* pj vidxulvi < p\dxul\oo\vi\'^ < p£^Cp\dxul\ 





2 




-dzVi 


=: Mie^ 


-dzVi 


£ 




£ 



Note that by Theorem 4.4, dxu\ G L°°{p). In the following, all the regularity results used in 
the estimates also follow from Theorem 4.4. 

★ For the next term, Poincare inequality is combined with Young inequality: 
P J V2dzulvi < p\dzul\oo \V2\ \vi\ < pCp\dzUl\oo \dzV2\ \dzVi\ 



<pc],\dzui\ 



-dzVi 

£ 



:M2 



* In a similar way: 

p u* ■ Vvi vi < pCp\ul\ 

J — 



-dzVi 



+ £^pCp\u*2\ 



:M3 



■■.Mi 



-dzVi 



Observe here that it was not possible to apply Lemma 5.1, since it* • n|an 7^ 0. 
• It remains the easier terms of Li and Ci (the ones which do not depend on v). 
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★ The first term is treated using again Poincare and Young inequalities: 
p / u*-Vu\vx < pCp\u*\^\Vu\\\dzVx\ < -{pCp\u*\oo\Vul\f+— 









-dzVi 




-dzVi 


e 




£ 



-dzVi 



★ Similarly, (1 - r)^ J d^ul vi<C+- 

★ The last term is estimated as follows, using Young inequality 

If 1 





2 




-OzVl 


<C + C 


-OzVl 


£ 




£ 



where c is a positive constant independent of e that can be chosen arbitrarily. 
Now, let us choose e and c small enough such that all constants satisfy: 

n M2£ Mge _ 9 (1 - r)v 

^'2'2' '2' - 36 



(5.5) 



Step 2. Let us multiply (5.1b) by £^V2 and integrate over J7. Again, the boundary terms in 
the integrations by parts vanish, since V2\dQ, = 0, and the convection terms are equal to zero since 
V • v = and v\Qa = (by Lemma 5.1). It follows: 



p^d_ 
2 dt 



\V2\'^ + il-r)iy\£VeV2\'^-^ J qdzV2 = -J Ti2dxV2-^J T22dzV2+ J L2V2 + ^J C2V2. 

Q. Q. n n n 



(5.6) 

Each term of / L2 V2 and / C2 V2 is estimated with the help of Poincare and Young inequalities 

n n 
as in the preceding step. 

2 



+ ^19^^21^ 1 + e'^ pCp\dzUl\oo \dzV2\''. 



■k £^p j V- Vu*2 V2 < £^ pCp I g^tt2 1 00 I \ 
Q =:M5 ^ 

•k £^p j U* ■ VV2V2 < £ pCp\ul\oo^ {\£dxV2\'^ + \dzV2\'^) + pCp\u2\oo |5zi^2p- 

fi =:Mr =-Ms 

* £^p J U* ■Vu*2V2<^£^p\u*\l,\Vu*2\^+£^^Cj>\dzV2\^ <C + £^Mg\dzV2\^. 



=:Mo 



★ By integration by parts (all boundary terms are equal to zero since V2\dn = 0) and Young 
inequality as before: 

{l-r)u£^J dlulv2 = -{l-r)u£^ j 5x«2 ^^^^2 < £ Q|9xW2l' + ^k5x^2p) ■ 
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* {l — r)iy / 9^n2 r2 < -^(1 — r)^z^^(7p|(9^tt2|^ + ci|52U2|^ < (7 + ci|92V2|^, where ci is a arbitrary 

J 4ci 
n 

positive constant. 

* / d^al2V2 < -^Id^c^ul'^ + ci\d^V2\'^ < C + ci\d,V2\'^ . 



* The C2 term is treated with integration by parts (again, no boundary terms since V2\dn 
vi\dn = 0) and the divergence equation. The term is then treated as the preceding one: 



^ J 9^(^22 V2 = ~ j 0-22 dzV2 = \ j ^22 ^xVl = ~~ j ^^'^22^1 



n n 

< Cp\dx(J22\ 



-OzVl 
£ 



n 

(J2 

< -r^\dxCT*22f + C2 

4C2 





2 




-dzVi 


< C + C2 


-OzVl 


£ 




£ 



Now, let us choose £, ci and C2 small enough such that 

-3 



|- , Mge^ , M7£, Mge, Mge, , ci , C2 < ■ 



(5.7) 



Step 3. After summing (5.4) and (5.6), and multiplying by 2rz/, it holds for e small enough 
(satisfying (5.5) and (5.7)): 

d 3 2ri^ f 

rvp-^^{\vi\^ + \ev2\^) + -r{l-ry {\VeVi\^ + \£VeV2\^)-^ / ^ (^^^i + ^^^2) < -V1-V2 + C, 

n 

where C is a constant independent of £. From the divergence equation V ■ v = dxVi + dzV2 = it 
follows that the pressure term J q {dxVi + dzV2) = 0, and equation (5.3) is obtained. 

□ 

Proposition 5.3. Let us suppose that 

y\dzUl\oo < 1/12, A*|cri2|oo < X, A*(|crii|oo + |<722|oo) < X, '^>^*\dz(^12\oo < X, A*|52(Tn|oo < X, 





1 


2 


1 


2 


1 






-m 


+ 2 


-T12 


+ 


-T22 






e 




e 




£ 





where x = — i^r(l — r). T/ien /or £ sma// enough, rn, ri2, r22 solution of (5.1) satisfy the 
following inequality: 

^^(|rnP + 2|n2|^ + |r22n+^ 

+ ^ (|V,riip + 2|Veri2p + |V,r22n < Pi + P2 + r(l - ry {\VeVi\^ + \£VeV2\^) + C, 

(5.8) 

where C is a constant independent of £. 

Proof. As in the preceding proposition, classical a priori estimates on rn, T12 and r22 are obtained, 
and the remaining terms are estimated accurately. 
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Til 

Step 1. Let us multiply (5. Id) by — and integrate over Q. Again, the convection terms 
/ V ■ Vrii Til contained in J dfTu th are equal to zero by Lemma 5.1, since V • v = and v\qci = 

(see (5.2)). Moreover, there is no boundary term in the integration by parts since the boundary 
conditions on cr have be chosen such that r • n|an = (see also (5.2)). It follows: 

2 



A* d 2 ^* 



J N{V,T12) 



Til + 



1 



-Til 

£ 

£ 



+ ^|V.riip 

e 



dxVlTll + - / LiiTll + ^ / -LiiTll. 

Q n 



• The terms of 



/a. 



Til are estimated as follows: 



n 





1 


' £2 


1 






-Til 


+ T 


-Tl2 


) 








£ 





::Mii 



\* f I £ 

★ In a same way: — / vidxCr\iTii < X*\dx(Jii\ooCp - 
e J V ' V ^ 

n =:Mi2 ^ 

Let us choose £ small enough such that: 

and 



-dzvi 



+ 



-Til 



Mii£^ < 1 

2 - 24 



^<Min|!ll^,l 
2 - 1 6 24 



(5.9) 



* A*y" 5a; ■U2 (712 Til < A*|cr^2|oo k5x'y2| 



1 

-Til 
£ 



1 

4C3' 



< A*|(Ti2|oo ( T— k^x'y2p + C3 



1 




-Til 




£ 





Here, it is not possible to choose C3 such that both coefficients are less than r(l — r)u/6 and 
1/24. Therefore, a condition on A*|(7j||2loo is imposed such that: 

< and A |ai2|ooC3 < -• 

Choosing C3 satisfying A*|c7i2|ooC3 = 1/24, the condition on A*|(Ti2|oo becomes: 



A*|o-i2|oo < --\/^(l - r) =: X- 



■k Similarly the following term can be estimated: 

V2dzCrliTii < A* 1^2(711 loo \dzV2\ 

CI 



£ 



-Til 



< >^*\dzCru\oo i ^1^2^21^ + C3 



1 


• 


-Til 




£ 





The same reasoning as before allows to control both terms providing that A*|5;2;(Tii|oo < X- 
. In order to treat the term -X* J u* • Vrnm, it is not possible to apply Lemma 5.1, since 
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u* ■ nloQ 7^ 0. However, integration by parts implies that 



-^7 



u* ■ Vrii Til 



u -riTii. 



On uj, since u* = (s,0) (see (4.2)), it holds u* ■ n = 0. Thus it remains to consider the 

boundary integral on F^. This boundary integral is split into two integrals on r+ and r_. On 

r_, it holds tt* • n > 0, thus — ^ / u* • nrf]^ < 0, and this term is trivially bounded by zero. 

r_ 

On r+, the boundary conditions are chosen in subsection 4.2 such that T|r+ = 0, therefore 

/ u* -riT? 

r+ 



rh = 0. 



• All other terms of 



Til are easier to manage, since they are linear in rn, and they are 



treated with Young and Poincare inequalities in a same way as the ones in vi, V2. 
For the terms of j L'n rn, we proceed as before: 



-Tl2 



-Til 



< X*\d,ul\ 



-Tl2 



+ 



-m 



Choosing A*|52;w5[||oo < 1/12, both terms are bounded by 1/24. 



A* f 
n 



-OzVl 
£ 



1 

-Til 

e 



'12|oo 



1 
4C3 



-dzVi 

£ 



+ C3 



1 

-Til 

£ 



Imposing A*|a"^2loo < X is enough to ensure that the coefficients are less than r(l — r)v/Q and 
1/24. 

2ti2 

Step 2. Now, multiplying equation (5.1e) by and integrating over fi, with the same 

reasoning as in the preceding step it follows: 



A* d I |0 A* f , , 

— 37 Tl2 -I / iV(v,rii - r22)Ti2 + 2 

£ at £ J 



1 

-Tl2 
£ 



+ ^|V.ri2|^ 

£ 



2rv 



= / dxV2 + -dzVl Tl2 + - / Li2 Tl2 + ^ / L'.o Tl2 



(5.10) 



The terms in L12 and L'12 are of the same type as the ones in Ln and L'^^^, and are treated very 
similarly to them, applying Young inequality, and assuming smallness assumptions on £. Thus, 
let us only write the terms needing additional assumptions. 



★ A* y d:cV2 (cTn - 0-22) T12 < A*(|crii|oo + |cr22|oo) \£dxV2\ 

a 

)<X. 
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-T12 



, and it is enough to assume that 



J 

J dzU* (rii - r22)ri2 < A*|5^'U*|oo ^ 



-Tl2 



* — / V2dzal2n2 <2X*\dzal2\oo\dzV2\ 
n 

A* 



that \*\dzUl\oc < 1/12. 
A* f 

* ^zVl{(Tu - 0-22) Ti2 < A*(|cr5Ll|oo + k22|oo) 



, and we assume that 2X*\dzal2\oo ^ X- 
1 





1 




1 


( 


-m 


+ 


-T22 




£ 







£ 



-T12 



1 

-T12 



it has already been assumed 



, it has already been assumed 



that X*{\al^\oo + Whloo) < X- 



T22 



Step 3. Multiplying (5. If) by — , and estimating the terms just as the ones in rn, it follows 



A* d , ,2 A* 
2£ dt £ 



J N(v, T12) T22 + 



-T22 



+ ^|V,r22p 

£ 



/ azW2 C + - / L22 T22 + ^ / -i^22 '^22- 



(5.11) 



Assuming that A|crj|"2|oo < X; A*|520"j|"i|oo < X and A*|52;'Ui|oo < 1/12, all the terms - / L22T22 and 



~2 / -^22 ''"22 are bounded and estimated as in Step 1. 



Step 4. Summing (5.9), (5.10) and (5.11), and noticing that 

-J N{v,Ti2)tii+ J N{v,Tii-T22)ti2+ J N{v,Ti2)t22 

= J (^£dxV2 - ■^^^^■"1^ (-'^12 ni + (rii - T22) ri2 + T12 r22) = 0, 



it follows that for £ small enough 



-Til 



+ 2 



-7-12 



+ 



-T22 



+ -(|V,Tii|2 + 2|Veri2p + |V, 



< Pi + P2 + r(l - ry {\VeVif + \eVeV2\^) + C, 
where we recognized the terms Pi + P2, and where C is a constant independent of e. 



□ 



Prom now on, let us come back to the notation with the superscripts denoting the dependence 
on £ and 77. 



Theorem 5.4. Suppose that the solution u*,a* of system (4.1)-(4.2) satisfies the following small- 
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ness assumptions 



y'idzulloo < 1/12, A*|c7i2|oo < X, A*(|crii|oo + |o-22|oo) < X, 2A* l^^a^s |oo < X, y\dzCru\oo < X, 

(5.12) 

where x = §y^r(l — r). Then the following convergences hold up to subsequences when rj and then 
£ tend to zero: 



Ui ul, dzUi dzUi, dxUi ^ dxU* in L {0,T,L (O)), 



0, dzU^^Q, dxU^^Q in L'^{0,T,L'^{n)), 



ecr 



^1 ^1) ^2 



in L'^{0,T,L'^{n)), 



0, £& 



ST] 



in L°°(o,r, L'^{n)). 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



Proof. Summing (5.3), (5.8) implies that for £ small enough (i.e. if assumption (5.12) is satisfied): 
d 



X* d 

:e 

£ 



erj 



+ \edxvir + \dzvT{']+^ 



eri\2 



1 



+ 



1 

-r- 



1 

+ 2 



1 

-r, 



< C. 



(5.17) 



Prom this inequahty, it follows that v'^'^' converges to v'^ in H^{Q) and r^'' converges in L^(r2), 
as Tj tends to zero, and are the solutions solutions of (5.1) without the terms tjAt^^. Indeed, 
recalling the weak formulation of the system (5.1), it suffices to notice that Holder's inequality 
allows to treat the term rjAr^'^: 

J V ' ^ 

n <c 

Moreover, and satisfy the following estimate: 



^dzvl 



e\2 







2 




2 1 








+ 




+ 2 


-r|2 



(5.18) 



It remains to pass to the limit as £ tends to zero. After integrating (5.18) between and T, it 

yields that 

lkillL2(L2) < ll^^if 111^2(^2) < Ce, thus the following convergences hold in Lp'{Q,T,L?'{p)) as 
£ tends to zero: 

^ and d^vl -> 0. (5.19) 
Prom these convergences, it follows that iif = + — > u\ in ^^(O, T, ^^(ri)) and dzu\ — > f^zuj 
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in L2(0,r, L2(n)). 

> ll^x^f ||L2(i2) < C, thus dxvf converges weakly in L^(0, T, Now, since it is already 
known that «f u^, it follows that dxuf dxU^ in L'^{0,T, L'^{^)). 

> Similarly ||f2llL2(L2) < l|f^2'i'2llz/2(L2) < C, thus ef^ and ec^zff converge strongly to zero in 

L2(0,r, ^2(1^)), and thus u| = e«| + £i;| ^ in L'^iO,T, L'^{Vl)), and 9^74 in L2(0,T, L2(J7)). 
C 

^ WxV2\\l'^{l'^) < — , thus dxU2 converges weakly in L^(0, T, L^(ri)). Since ^1 — >■ 0, it implies 
that dxul ^ in L'^{0,T,L'^{n)). 

» lkiillL2(L2), lki2llL2(L2), ||r|2||L2(i2) < Ce, therefore rfi,rf2,r|2 ^ in ^^(0, T, L2(f2)). 
Thus eafi = o-^i+rf^ ^ o-jfi in L'^{0, T, L'^{^)), and in the same way eal2 (^12 ™ ^! ^"^i^))^ 
eal2^a*22 in L2(0, T, L2(Q)). 

> Prom the terms with the derivatives in time, using the fact that v^\t=o = Uq — u* E -L^(17) 
and T^\t=o = 0*0 — (7* G are bounded independently of e, we can conclude that 

II'*^^IIl°°(l2) < C and ||'''^||L°°(i2) < C^/e. 

These estimates and the uniqueness of the limit imply that vf and eff converge weakly-* in 
L°°(0, T, toward zero, and that converges strongly in L°°(0, T, L^(ri)) toward zero, 

which proves the last estimate (5.16). 

□ 

Note that in a simplified case (with a simpler geometry), the hypothesis (5.12) is satisfied under 
a small data assumption on the physical parameters. 

Remark 5.5. When h is constant with respect to x, p* is also independent of x, so that equation 
(4.1) reduces to 

-(i-r)a.v-4 i . , '::L..^ =^- 



dz Vl + A*2|5^n*|2^ 

It has been shown in [8] for example that for r < 8/9 this equation admits a unique solution 
«t = .(l-f). 

Now, it follows that a*2 = i ^ x*'^\dlu\\'^ ^ h + X*^s^/h ' ""'^ ^^^^ ^ ""^^2 = -A*5^M*a*2 = 
-rus'^X* 

In this case, hypothesis (5.12) becomes more simple. Since dzu\ = —s/h, a\-y and 0^2 o,re constant 
with respect to z, so that the last two conditions are trivially verified. Using the fact that r < 8/9, 
it leads to a smallness condition on sX* with respect to h (sX* < h/12 is enough in order to satisfy 
all conditions). 

Observe that this condition is not optimal, but it shows that in the simplified case when h{x) is 
constant, a simple choice of the parameters s, X* and h satisfies hypothesis (5.12). 

5.3 Convergence of the pressure 

It remains to prove the convergence of the pressure. 
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Theorem 5.6. Under the same smallness assumption (5.12), the following convergence result 
holds up to a subsequence for p: 



Proof. Throughout the proof, C will denote some generic 



■^p ^ p* inV'{0,T,L'^{n)). (5.20) 

,iie pruui, will ueiiuue some geiieni; constants independent of e. Let £ < 1. 
Let us integrate over J7t = x (0,T) equation (5.1a) multiplied by e'^ipi, for any function 01 G 
H^{n). It follows: 



r^T^ 






(5.21) 



Using the fact that (pi is independent of t, the first term becomes 

T 

pe^y dtvict>i = pe^ j <Pi J dtvi = pe^ J Mn{T) - viiO)), 
Qt no n 

where vi{0) = — -u^ denotes the value of vi at time t = 0. Now, introducing 

T 

TT = y Q'dt, 



and using integration by parts for the pressure term (the boundary term is zero since (pi G Hq{0,)), 
(5.21) becomes: ^cpi G Hq{Q), 

pe^ J (f)i{vi{T) - vi{0)) + pe^ J vidxVi(f)i + J r2dzVi(f>i + (1 - r)u£^ J da:Vi d^cpi 

: J Tiidx(l)i- J Ti2dz(t>i + e^ J Li^i + e J Cicpi. 
fijH ri'j-' fi'j^ ri'j^ 



+ 



(1 - ry j dzVi dzcpi - j -Kd^cpi = -e 



It remains to estimate all terms independent of tt. The non-linear terms are to bee handled with 
care, since <pi ^ Lr^iyi). Proceeding as in [3], Holder inequality with exponents 2 + 5, (5' and 2 
leads: 



y VldxVl(pl <\(f)l\s' J \vi\2+S \dxVl\, 



(5.22) 
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1 11 ^ 2(2 + (5), , 
where - — - ~^ 2 ~^ ~ ^ (which implies that 6 = ). According to interpolation theory, 

[L^jL^lfl = L^"'"'^ for 6 = -, and the following estimate holds: 

L le 2 + S 



\vi\2+s < C\vi\l\vi\^-'^ . 
Moreover Lemma 3.2 of [1] states that for vi G Hq{Q,), it holds: 

\viU<V2\d,vi\'/''\d,vif/\ 
Using the two last inequalities and Poincare inequality, (5.22) becomes 



pe 



Q7 



and Holder inequality implies that 

pe^ J vid^vi(f)i < pe'^\4>i\s'C\\da;Vi\\\t^^l^l) 11^2^111^2 



i-e/4 



Now, choose 6 (and thus 5) such that S' > 6. It suffices to take < |, for example take 6 = -. 
Then S' = 6, and the usual Sobolev embeddings read H^{Q,) L^(ri) (which is true in dimension 
2 or 3). Therefore, the last estimate becomes 



pe'^ J vid^vi(t)i < pe^C\\(t)i\\Hi\\dxVi\\fJ^^^-^ 



,,11/12 



Now, recalling that H^^?;! 11^2(^2) < Ce and Ht^ifi || l2(l2) < C, we conclude 

pe^ J vAvkPi < pe^CUiWHie^'/'^ = pe''+'^/^''C\\<Pi\\Hi < Ce\\Mm- 



In a similar way, it holds 



pe J r2d,vicj)i < ps'^-^^^^CUiWh^ < Ce||</>i||^i. 



For the term ps^ J (pi{vi{T) —vi{0)), we apply Cauchy-Schwarz inequality. t;i(0) is bounded, and 

a, 

for vi(T), we use Poincare inequality. It follows, using the fact that l^^-uil < Ce: 

pe'^J MviiT)-v,{0)) < iC\vi\+C)e''\\MH^ < iC\d,vi\+C)e^\\MH^ < Ce^cPiWjj, < C£||<^i||^^i. 
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For the other hnear terms, a simple apphcation of Cauchy-Schwarz inequahty allows to obtain 
similar estimates. Indeed, it suffices to use the estimate (5.18) in order to estimate the L^-norm 
oi dxVi, dzVi, Til, ri2, Li, Ci. For example, since \dxVi\ < C, the following estimate holds: 

pe"^ J dxVidx(j)i < pe'^\dxVi\ \dx4>i\ < C£^||^i||iji. 

For the terms Li and Ci, Ci and the constant part of Li are obviously bounded uniformly in 
£. It remains to estimate the linear term Ci of Li. Recalling its definition and using Poincare 
inequality in the second estimate: 

\Ci\ <C{\vi\ + \V2\ + \dxvi\ + \d,vi\) <Ci\d,vi\ + \dxvi\ + \d,V2\) . 

Using again (5.18), the boundedness of £1 follows: 

Hence ^cpi G H^in): 

J dxTr(f)i <C {e + e^\dxvi\ + + £|rii| + |ti2| + + £|Ci|) ||(^i||^i < C£\\(f)i\\Hi. 

The same approach with (5.1b) gives a similar estimate, for all (1)2 G Hq{Q): 

J dzTr(p2 < C [e + e'^\dxV2\ + e^\dzV2\ +£^|ti2| +£|t22| +£^|i^2| + £|C2|) WhWn^ < Ce\\(j)2\\Hi- 
n 

Thus we can conclude that || V7r||/,oo(j:/-i) < Cs. 

Now recall that for / G -^^o(^)' holds that |/| < ||Vg||j:/-i (see for example [17]). Since p G I/o(^) 
and p* G Lq (J7) , q lies in Lq (J2) . From the definition of tt as function oiq, it is clear that tt G Lq (J7) . 
This allows to deduce 



:n L-^[0, '1\L^q[\1)) when £ ^ U. iNow, since q= ^ 
to zero in P'(0,r,L^(O)), and therefore: 



thus TT tends to zero in L°°{Q. T,L'^{^)) when £ ^ 0. Now, since q = — , it follows that q tends 



e'^p ^ p* in P'(0,T,l2(O)). 

£— >0 

This finishes the proof. □ 
5.4 Open problems 

This work concerns only the solutions of the problem (3.1) that are obtained as the limit of the 
regularized problem we chose (with an additional term — ryAtr). Since there is no uniqueness result 
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for problem (3.1), it is not known how other solutions behave. 

Formally, the passing to the limit can be done for a 7^ (see [4]), and a similar limit problem 
(involving the parameter a, but of the same structure). However, the proof of the existence 
theorem in Cl^ strongly relies on the fact that o = 0. No global results are proved in the case 

Last, since the computations are independent of the dimension of the domain Q, the result should 
be true in the three-dimensional case. The limit problem on (tt*,^*,^*) reads: 

' (1 - r)i^d^,ul - d,p* + 5,a*3 = 0, 
(1 - r)vdlu2 - dxP* + dza2s = 0, 

d,p* = 0, 

V • u* = 0, 

-Xy.ulat^i + a*u = 0, 

-A* 9^^2(723 + 0-22 = 0, 

A* A* 

Y^zuKal^ - 0-22) - Y^^'^*^'^^'^ + ^23 = rudzu*2, 

A* {dzu\al^ + 9^^20-23) + 0-33 = 0, 
A* A* 

y 9^Mi(c733 - c7n) - Y5^W2f^^2 + '^12 = rudzu\. 
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